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Abstract. The classical hard core model from statistical physics, with 
activity A > and capacity C = 1, on a graph G, concerns a probability 
measure on the set 1(G) of independent sets of G, with the measure 

F^""* . of each independent set I £ 1(G) being proportional to A' 7 '. Ramanan 

Cn ■ et al. proposed a generalization of the hard core model as an idealized 

model of multicasting in communication networks. In this generaliza- 
tion, the multi-state hard core model, the capacity C is allowed to be 
a positive integer, and a configuration in the model is an assignment of 
states from {0, . . . , C} to V(G) (the set of nodes of G) subject to the 

, <""| constraint that the states of adjacent nodes may not sum to more than 

C. The activity associated to state i is A", so that the probability of a 
configuration a : V(G) -»• {0, . . . , C} is proportional to X^^v(G) ct M. 

In this work, we consider this generalization when G is an infinite 
rooted &-ary tree and prove rigorously some of the conjectures made by 
Ramanan et al. In particular, we show that the C — 2 model exhibits a 
^ ' (first-order) phase transition at a larger value of A than the C = 1 model 

\& . exhibits its (second-order) phase transition. In addition, for large b we 

^^ ■ identify a short interval of values for A above which the model exhibits 

phase co-existence and below which there is phase uniqueness. For odd 
C, this transition occurs in the region of A = (e/6) 1 ' ' ' 2 ' , while for even 

C*"*- ■ C, it occurs around A = (log b/b(C + 2)) 2/(c+2) . In the latter case, the 

^^ ' transition is first-order. 
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1. Introduction 

1.1. The Multi-State Hard Core Model. Let G = (V, E) be a finite or 
countably infinite graph without loops, and let S be a finite set. We refer to 
the elements of S as states. Many stochastic processes on S that arise in 
applications are subject to "hard constraints" that prohibit certain values of 
S from being adjacent to one another in the graph G. Such processes only 
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attain configurations that lie in a certain feasible subset of S . A generic 
example is the hard core model, which has state space S = {0, 1} and im- 
poses the constraint that no two adjacent vertices in the graph can both have 
the state 1. In other words, the set of feasible configurations for the hard 
core model on a graph G is {a € {0, 1} V : a x + o~ y < 1 for every xy G E}, 
or, equivalently, the collection of independent sets of the graph G. Pro- 
cesses with such hard constraints arise in fields as diverse as combinatorics, 
statistical mechanics and telecommunications. In particular, the hard core 
model arises in the study of random independent sets of a graph [5, 7], the 
study of gas molecules on a lattice [2] , and in the analysis of multicasting in 
telecommunication networks [9, 11, 15]. 

In this work, we consider a generalization of the hard core model, which 
we refer to as the multi-state hard core model, in which the state space is 

S c = {0,1,2,..., C}, 

for some integer C > 1, and the set of allowable configurations is given by 

^G = {c G Sq '■ cr x + a y < C for every xy € E}. 

When G is the d-dimensional lattice TL d , this model was introduced and 
studied by Mazel and Suhov in [14], motivated by applications in statistical 
physics. In our work, we focus on the case where G is an infinite rooted 
6-ary tree (i.e., an infinite graph without cycles in which each vertex has 
exactly 6+1 edges incident to it, except for one distinguished vertex called 
the root which has b edges incident to it), which we denote by T b . 

On the tree, this model was studied by Ramanan et al. in [15] as an ideal- 
ized example of multicasting on a regular tree network, each of whose edges 
has the same capacity C. In communications, multicasting arises when, in- 
stead of having a simple end-to-end connection, a transmission is made from 
a single site to a group of individuals [1]. An important performance mea- 
sure of interest is the probability of packet loss for a given routing protocol 
[17]. As in [15], here we consider an idealized model in which the routing is 
simple in the sense that nodes multicast only to their nearest neighbors, and 
study the impact of the connectivity of the network (i.e., the value of 6) and 
the arrival rate on the blocking (or packet loss) probabilities. The state a v 
of any node or vertex v £ V represents the number of active multicast calls 
present at that node. Multicast calls are assumed to arrive at each node 
as a Poisson process with rate A and require one unit of capacity on each 
of the 6+1 edges emanating from that node. If this capacity is available, 
then the call is accepted and the number of active multicast calls at that 
node increases by one, while if the required capacity is not available, then 
the state of the node remains unchanged and the call is said to be blocked 
or lost. Calls that are accepted require a random amount of service and 
then depart the system. Service requirements of calls are assumed to be in- 
dependent and identically distributed (without loss of generality with mean 
1), and independent of the arrival process. This model is a special case of 
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a loss network (see [9] for a general survey of loss networks and [12, 15] for 
connections with this particular model). 

For a finite graph G and arrival rate A, it is well-known that the associated 
stochastic process has a unique stationary distribution hg,\ ° n ^G that is 
given explicitly by 

(1) MG,A = ~ — ]T A CT „ for a G Q G , 

where Zg,\ is the corresponding normalizing constant (partition function) 
%G,\ = J2aen G YlveV ^v > where the form of Aj depends on how the multicast 
calls are served. If the calls are assumed to be served in a first-come first- 
served manner at each node (see [9]), then we have 

A* 
Aj = — , i = 0,...,C, 

i\ 

If they are served using the processor sharing scheduling discipline at each 
node (see [10]), then we have 

(2) \ = \\ i = 0,...,C. 

Here we adopt the usual convention that 0! = 1, so that Ao = 1 in both 
models, and we will sometimes refer to the arrival rate A as the activity. In 
this paper (as in [14]) our Aj's will always be as defined in (2). Thus, our 
exclusive focus will be the study of the multi-state hard core model on a 
6-ary tree T b with activities given by (2). 

1.2. Gibbs Measures and Phase Transitions. Although there is an ex- 
plicit expression (1) for the stationary distribution on a finite graph, the 
computational complexity of calculating the normalization constant for large 
graphs limits the applicability of this formula. Thus, in order to gain insight 
into the behavior of these measures on large graphs, it is often useful to con- 
sider the associated Gibbs measure on an infinite graph. Roughly speaking, 
a Gibbs measure on an infinite graph G associated with an activity A is char- 
acterized by the property that the distribution of the configuration on any 
finite subset U of V, conditioned on the complement, is equal to the regular 
conditional probability of the measure £*G[0i a on the restriction G[U] of the 
graph G to the closure U = U U dU of U, given the configuration on the 
boundary dU of U (see Definition 2.1 below for a more precise formulation). 
It is not hard to show that such a Gibbs measure always exists (in a far 
more general context, see for example [8]). 

However, unlike stationary distributions on finite graphs, the associated 
Gibbs measures on infinite graphs may not be unique. If there are multiple 
Gibbs measures associated with a given arrival rate or activity A, we say 
that there is phase coexistence at that A. Let T n denote the finite sub-tree 
of T with root r and depth n, which contains all vertices in T that are at 
a distance of at most n from the root r. As is well known (see, for example, 
Chapter 4 of [8]), for a fixed activity A > 0, one way to obtain a Gibbs 
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measure on the tree T rooted at r is as the suitable limit of a sequence of 
measures, where the nth measure in the sequence is the stationary measure 
UT n uc9T n ,A on T n U dT n (as defined in (1)), conditioned on the boundary dT n 
being empty (i.e., conditioned on all vertices in the boundary having state 
0). We shall refer to this Gibbs measure as the empty boundary condition 
(b.c.) Gibbs measure (corresponding to the activity A). In a similar fashion, 
we define the full b.c. Gibbs measure to be the limit of a sequence of con- 
ditioned measures on T n , but now conditioned on the boundary dT n being 
full (i.e., conditioned on all vertices in the boundary having state C). Let S\ 
denote the total variation distance of the marginal distributions at the root 
r under the empty b.c. and full b.c. Gibbs measures corresponding to the 
activity A. When A lies in the region of uniqueness, clearly the empty b.c. 
Gibbs measure coincides with the full b.c. Gibbs measure, and so S\ = 0. 
On the other hand, when A is in a region of phase coexistence, then 5\ > 
and it can be shown (due to a certain monotonicity property of our model 
established in Lemma 2.2 and Proposition 2.3) that the empty b.c. and full 
b.c. Gibbs measures must necessarily differ. If there exists A cr = X cr (C) for 
which there is uniqueness for each A < A cr and phase coexistence for every 
A > A cr , then we say that a phase transition occurs at A cr . Moreover, if 
5\, as a function of A, is continuous at A cr , then we say that a second- order 
phase transition occurs, while if 5\ is discontinuous at A cr , then we say that 
a first-order phase transition occurs. 

When C = 1, the phase transition point A cr (l) on the tree is explicitly 
computable and is easily seen to be a second-order phase transition (see 
[9, 18, 19] and also Section 2.2). The behavior is more complicated for 
higher C. The multi-state hard core model on the d-dimensional lattice 
Z was studied in [14], where it was shown that when C is odd, there is 
phase coexistence for all sufficiently large A, while when C is even, there is a 
unique Gibbs measure for each sufficiently large A. If phase coexistence were 
known to be monotone in the activity (this remains an open problem on Z 
even when d = 2), then the result of Mazel and Suhov would imply that 
there is no phase transition on Z for even C . On the other hand, numerical 
experiments for the multi-state hard core model on the regular tree (see 
Section 3.5 and Figure 5 of [15]) suggest that there is a phase transition on 
the tree for every C, but that the order of the phase transition depends on 
the parity of C (being first-order for even C and second-order for odd C). 
This is particularly interesting as it shows that the parity of the capacity has 
an effect on the regular tree as well, although the effect is not as pronounced 
as on the d-dimensional lattice. 

The study of phase transitions of models with hard constraints on trees 
has been the subject of much recent research (see [4, 5, 13]). In [4], the 
focus is on classifying types of hard constraints (as encoded in a so-called 
constraint graph) on the basis of whether or not there exists a unique sim- 
ple invariant Gibbs measure for all activity vectors (\i,i £ S). For C > 1, 
the model that we present here allows for two l's to be adjacent, but never 
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allows a 1 to be adjacent to C which, in the language of [4], implies that 
the associated constraint graph is fertile. From Theorem 8.1 of [4] it follows 
that there exist some activity vectors for which there exist multiple simple 
invariant Gibbs measures. However, the emphasis of our work is quite dif- 
ferent, as our aim is to identify regions where multiple Gibbs measures (not 
necessarily simple and invariant) exist for the particular choice of activity 
vector given in (2). Another related work, again motivated by telecommu- 
nication networks, is [13], which studies Gibbs measures associated with a 
three-state generalization of the hard core model. However, the hard con- 
straints considered in [13] are somewhat different from the C = 2 case in 
our model. 

1.3. Main Results and Outline. The main contribution of this paper is 
to make rigorous some of the conjectures made in [15], leading to a bet- 
ter understanding of the multi-state hard core model. Our results may be 
broadly summarized as follows. 

(1) For C = 2 and every b € N, b > 2, we show that the Gibbs measure 
is unique for larger values of A than in the usual C = 1 hard core 
model (see Corollary 3.2) and we also show that the phase transition 
is first-order (see Theorem 3.3). Recall that, in contrast, for C = 1, 
the phase transition is second-order. 

(2) For large values of b, we identify a rather narrow range of values for 
A, above which there is phase co-existence and below which there is 
uniqueness. Although we do not establish the existence of a unique 
critical value X cr (C) at which phase transition occurs, we establish a 
fairly precise estimate of A cr (C) if (as we strongly believe) it exists: 
when C is odd, 

(see Theorem 4.1), while for C even, 

1 log b 



A C r(C) 



C + 2 b 



2 

C + 2 



(see Theorem 4.5). 
(3) For all even C and all sufficiently large b (depending on C), the model 
always exhibits a first-order phase transition (see Section 4.3). 

The outline of the paper is as follows. First, in Section 2 we establish a 
connection between phase coexistence and multiplicity of the fixed points of 
an associated recursion. This is based on the construction of Gibbs measures 
as limits of conditional measures on finite trees with boundary conditions, as 
mentioned above. In Section 3 we provide a detailed analysis of the recursion 
in the special case C = 2. In Section 4 we study the recursion when b is 
large and identify the phase transition window. Finally, in Section 4.3 we 
study the asymptotics for large b when C is even and provide evidence of 
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a first-order phase transition. An interesting open problem is to rigorously 
establish that the phase transition is second-order for all odd C. 

2. Gibbs Measures and Recursions 

2.1. Gibbs Measures on Trees. Consider any graph G = (V,E) with 
vertex set V and edge set E C V^ 2 ' (the set of unordered pairs from V). For 
any U CV, the boundary of U is dU = {x G V\U : xz G E for some z G U} 
and the closure of U is U = U U dU. Let G[U] denote the restriction of the 
graph to the vertex set U. For a G S 1 ^, let ay = (a v ,v G [/) represent the 
projection of the configuration a onto the vertex set U. With some abuse 
of notation, for conciseness, we will write just a v for at v \ and refer to it as 
the state or, inspired by models in statistical mechanics, the spin value at 
v. For U C V, let J~(U) be the <7-field in Sq generated by sets of the form 
{a v = i} for some v G U and i G Sc- We now provide a rigorous definition 
of the Gibbs measure. 

Definition 2.1. A Gibbs measure for the multi-state hard core model as- 
sociated with the activity X is a probability measure [i on (Sq, J-(V)) that 
satisfies for all U C V and fi-a.a. r G S)j, 

fl(ay = Ty\a v \y = Ty\y) = ^ G [U],x( a U = Ty\a d y = Tgy), 

where fJ>G\u] X ^ s as defined in (1), with \ given as in (2). 

We now specialize to the case when G is a regular, 6-ary, rooted tree T b 
with root r. A child of a vertex x in T b is a neighboring vertex that is further 
from r than x; the vertices (other than x) that lie along the path from x to r 
are the ancestors of x. We will be concerned with (complete) finite sub-trees 
T of T b rooted at r; such a tree T is determined by a depth n, and consists of 
all those vertices at distance at most n from r. It has \T\ = (b n+l — 1)/(6— 1) 
vertices, and its boundary dT consists of the children (in T b ) of its leaves 
(so that \&T\ = b n+1 ). The tree consisting of all vertices at distance at most 
n from the root r will be denoted by T n . 

Given a finite sub-tree T and r G i7 T b, we let QTp denote the (finite) 
set of spin configurations a G ^tu9T that agree with r on dT; thus r 
specifies a boundary condition on T. For a function / : QtuOT -> twe 
denote by A*ta(/) = So-en r /'taW/W ^ ne expectation of / with respect 
to the distribution /ij, A (<r) ex Y\ veT X av . On the configuration space Qjb, 
we define the partial order a -< n if and only if a v < r\ v for all v with even 
distance d(v,r) from the root and a v > rj v for all v with odd distance from 
the root. Given two probability measures on f2 T b, we then say that /x -< v 
if //(/) < v(f) for any (bounded) function / that is non-decreasing with 
respect to the above partial order. 

Let T be a complete finite tree rooted at r, and let /ij, A and /j,j, A , respec- 
tively, be the empty b.c. and full b.c. measures (corresponding to the two 
boundary conditions identically equal to and C, respectively, on dT). The 
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following monotonicity result is well known (see, for example, Theorem 4.1 
of [19]). However, for completeness, we provide an independent proof of this 
result, which involves a Markov chain argument that constructs a simul- 
taneous coupling of (Mt'A'^ta'Mt'a) such that the required monotonicity 
conditions are satisfied with probability one. 

Lemma 2.2. For any r G Qjb, 

lAr x "^ Mt a "^ Mr a if d(dT,r) is even, 
Hj,\ ~< Mr a "^ Mr a if d(dT,r) is odd. 

Moreover, if d(dT, r) is even (respectively, odd) there is a coupling ttt = 
(a°,a T ,a c ) of (/Jj, x , /i T A , /j,j, ^) such that a -< a T -< a c (respectively, a c -< 
a T ~< a ) with probability one. 

Proof. We consider only the case when d(dT, r) is even, since the other case 
can be established in an exactly analogous fashion. On O^ x fij, x Qj, we 
construct an ergodic Markov chain {a°(t), cr T (t), o - (£)}t e z + such that at any 
time t € Z + the required ordering relation cr (t) -< o~ T (t) -< a c '(£) is satisfied, 
and moreover each replica is itself an ergodic chain that is reversible with 
respect to the measure [i T x with the corresponding boundary condition. 
The stationary distribution ttt of the global chain will then represent the 
sought coupling of the three measures. 

The chain, a standard Heat Bath sampler, is defined as follows. Assume 
that the three current configurations corresponding to 0, r and C boundary 
conditions are equal to (a, (3, 7) respectively and that they satisfy the order- 
ing relation. Pick uniformly at random v £ T and let (a, b, c) be the maxi- 
mum spin values in v compatible with the values of (a, /?, 7) on the neighbors 
of v respectively. Due to the ordering assumption either c < b < a or the op- 
posite inequalities hold. Then the current three values at v are replaced by 
new ones, (a' v , /3' v , j' v ), sampled from a coupling of the three distributions on 
{0, 1, . . . , a}, {0, 1, . . . , b}, {0, 1, . . . , c} which assign a weight proportional 
to A* to the value i. It is clear that such a coupling can be constructed in 
such a way that (a' v , /3' v , 7^,) satisfy the opposite ordering of (a, b, c) and thus 
the global ordering is preserved. □ 

Consider now the sequence {T2 n } n eN with d(dT2 n ,r) = 2n. Then, thanks 
to monotonicity, linin^oo fij, \ = fJ-\ exists (weakly) and it defines the 
maximal Gibbs measure. Similarly linin^oo ^ a = Ma defines the minimal 
Gibbs measure [8]. Notice that, by construction, lim n _ s>00 fij, a = Ma 
while linin^oo \£j 2n \ = H\- Finally, for any other Gibbs measure fi, it 
holds that /i° -< n\ ■< jjf^ . 

The main problem is therefore that of deciding when ^ = Mx- ■'■ n what 
follows we establish the following criterion, which is in fact an equivalent 
criterion, since the other implication is obviously true (see [18, 19]; see also 
[3, 16] for a similar discussion in the special case of C = 1). Let P^ A be the 
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marginal of (ij, A on a r given boundary condition r, and let P A and P A be 
the corresponding marginals for /j,*? and //9, respectively. 

Proposition 2.3. For every A > 0, ifP% = P A then /j,*? = /j,^ . 

Proof. Assume P A = P A . Then, by monotonicity, 

(3) lim ||P£ A -P° jA || TF = 0, 

where || • \\tv denotes total variation distance. Let A be a local event 
(i.e., depending only on finitely many spins) and let m be sufficiently large 
so that A does not depend on the spin configuration outside T m . Fix n > m, 
and let 7T2 ra = (o-°,a T ,a c ) be the monotone coupling of (/j, j, A ,/U^ A'^t a) 
described in Lemma 2.2. Then 

\^T 2n ,x( A ) ~ /"L,a(^)II < ^(o). T a\>. for souic r G T, 

o 



'r 2n ,\( A ) ~ VT 2n , 

c 



^ E X>2n(<tf > fc ><r2) 

ver m k=0 

d(v,r) even 

c 

+ E E^ 2 "^ > k > a°) 

vGT m k=0 
d(v,r) odd 

c 

E E [»T 2n ,\(<7v >k)- »T 2n , X (o- v > k)] 

v£T m k = 

d(v,r) even 

c 

d(v,7") odd 

For simplicity, let us examine an even term /j,j, x (a v > k) — fij, x (a v > k) 
and show that it tends to zero asrn oo. Let w be the immediate ancestor 
of v. By conditioning on the spin value at w we can write 

A*t 2 „,aK = *) - Mr 2n ,AK = *) 

C-i 

= E ^T 2n ,\{^w = j) [mt 2 „,aK = « k«» = J) - mt 2 „,aK = * kt» = ;')] 

j=0 

C-i 

+ E [^T 2 „,a(^ = J) ~ Mt 2 „,a(^ = J')] Mt 2 „,aK = i\cr w =j). 
3=0 

By iterating upwards until we reach the root, and using (3), we see that it 
is enough to show that 

lim max max max \u Tn \{o~ v = i I o~ m = j) — u Tn x (<r„ = i I o~ w = j)\ = 0. 
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Now, let 

(4) 



Z T k ii) 



x 1 yi n xav 

CTef2 T fe \{r}^T fe \W 



denote the partition function (or normalizing constant) on the complete 
finite tree Tk with boundary conditions r and a r = i. It is clear that 



zm 



Therefore, 



^T 2 „,a(^ = i\cr w =j) 



Pfc,A(0 

n,A(o) ^ r (o) ' 

Z^k<C-j ^2n-n„ W 



2n-n„ W 



^n-nJ^C-j) 



where n^ denotes the level of v (the distance from the root). A similar 
relation holds for the full boundary condition. 

The proof is concluded once we observe that n v < m and that 



3>0 

In— n 



v ,x(<?r<C-j)>Fl x (0)>0. 



D 



2.2. Recursions. Our next step, as in many other spin models on trees, 
is to set up a recursive scheme to compute the relevant marginals P° » and 
P„ ^. In what follows, for simplicity we count the levels bottom-up and the 
boundary conditions are at level 0. Moreover, since the recursive scheme is 
independent of the boundary conditions, and since we will never be consid- 
ering more than one value of A at a time, we drop both from our notation. 
For i = 0, . . . , C, and n £ N, we set 



Qn{i) ■-- 



Pn(0) 



l \ Mi) ■= S-°to W( ,!i = I 1 ~ P (<* > C ~ i)]- 1 . 



Z^k=0 Qn(fc) 



Thus Rn(0) = 1 and R n (i) < R n (i+1). Moreover, let Z n be as defined in (4), 
but with r equal to the empty b.c.Then we obtain the recursive equations 



(5) 



Z n +l(i) 



Qn+l(i) 



Rn+l(i) 



E 

E 



C-i 

J2Zn(k) 
,fe=0 

Efc=Q Qnjk) 

z^k=oQ n (k) 

• \ k 
fc=0 7^(fc) 



RiiiY 



c- 



\ k 



fc=0 Rb n (k) 
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The case when C = 1 (the usual hard core model) can therefore be studied 
by analyzing a one-dimensional recursion governed by the following maps: 



(6) Ax) ■= j^—T- h , J 2 (x) := J(J(x)) 



(1 + x) 6 ' ^ ' v y " (1 + 



(l+x) b 



Indeed, J defines the recursion for the quantity Z n (l)/Z n (0), while Ji defines 
the recursion of this quantity between two levels on the tree. We close this 
section with a summary of the properties of J and J 2 which, when combined 
with Proposition 2.3, show that A cr (l) := b b /(b— l) b+l is the phase transition 
point for the standard hard core model (see, for example, [9]), and that the 
phase transition for C = 1 is second-order. These properties will turn out to 
also be useful for our analysis of the higher-dimensional recursions {i.e., when 
C > 2). We start with the definition of an S-shaped function. 

Definition 2.4. A twice continuously differentiable function f : [0, 00) 1— >■ 
[0, 00) is said to be S'-shaped if it has the following properties: 

(1) it is increasing on [0, 00) with /(0) > and sup x f{x) < 00; 

(2) there exists x € (0, 00) such that the derivative f is monotone in- 
creasing in the interval (0,x) and monotone decreasing in the inter- 
val (x, 00); in other words, x satisfies f"(x) = and is the unique 
inflection point of f . 

For future purpose, we observe here that the definition immediately im- 
plies that for any 8 > 0, and S'-shaped function /, Of is also an S-shaped 
function. It is also easy to verify that any S-shaped function has at most 
three fixed points in [0, 00), i.e., points x € (0, 00) such that f{x) = x. We 
now summarize the salient properties of J2 (see e.g. Fig. 1), all of which 
may easily be verified with some calculus. 

(1) J2 is an S-shaped function with J2(0) = A/(l + \) b and sup^, J2{x) = 
A, and a unique point of inflection x* € (0, 00). 

(2) J has a unique fixed point, xq, which is also a fixed point of J2. 

(3) If A < A cr (l) then J^ix) < 1 for any x > and xq is the unique fixed 
point of J2. 

(4) If A > A cr (l), then J2 has three fixed points x_ < xq < x + , where 
J(x-) = x + and J{x+) = x_. Moreover J 2 {xq) > 1, J^ix) < 1 
for x G [0,x_] U [x + ,+oo) and the three fixed points converge to 
zo(Acr(l)) as A I A cr (l). 
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The supercritical function J_2 




FIGURE 1. Graph of the function J2&) for b = 2, A = 7 
(A cr = 4) 



3. Analysis of the recursions when C = 2 
When C = 2 we have i? n (l) = [l — ¥ n (a r = 2J] and (5) can be written 



as: 



(7) 



Rn(0) = 1, 

R n+1 {2) = 1 + ^^+ Xr 



Rn+iO-) 



Ri(i) Ripy 



1 + 



Mm 



1 + 



WX) 



On replacing n by n — 1 in the last equation above, we see that 

A \ 



Rn(2) = R n {l) 1 + 



R-n-ii 1 ) 



Substituting this back into (7), we obtain an exact two-step recursion for 
Y n := R n {l): 

(8) Y n+1 = 1 + 




[Yb + \] 
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It is useful to determine the initial conditions (Yq, Y\) for the recursion given 
the boundary conditions at the level. 



Oo^i 



(+oo,l) 



if the b.c. is full (i.e., identically C) 
if the b.c. is empty (i.e., identically 0) . 



Numerical calculations of (8) using Mathematica strongly suggest that the 
critical value X cr , below which the recursion settles to a limit independent 
of the initial values, takes approximately the following values: 



b 


A cr 


2 


7.2753875 


3 


3.58029 


10 


1.107665 


100 


0.2817409 



and that the transition is always first order (i.e., if limsup n Y n 7^ liminf n Y n 
then their difference is strictly larger than some positive constant 5). Similar 
observations were made in [15] (see Section 3.4 therein). Here, we provide 
a rigorous proof of these results. 

Let us change variables from Y n to X n := Y n — 1 in (8). It then follows 
that 

A 2 In) 

••;■+-! ^ 7 77 ZJ = F+ {X n -l), 



(9) 
(10) 



X n 



min J>n Y j + ^ 



1 + 



A 



(1+X n ^) b 



X 



n+l 



> 



maxj>„ Yf + A 



1 + 



A 



(1+X„_i)* 



F^(X n ^) 



M 



The maps F± defined above can be rewritten in terms of the map J2 defined 
in (6) as follows: 



W, 



F 



F[ n) (x) 



maxj> n YP + X 



min i>n Y j + ^ 



J 2 (x). 



A 



-Mx) 



Next, for k > 0, we define 

(11) F.(,):. k + a . 

so that Fq = J2. For any k > 0, F K is a strictly positive multiple of J 2 
and hence is also an S'-shaped function (with the same inflection point x*). 

If we denote the fixed points of F K by x_ ' < Xq < x^' (with the obvious 
meaning) we see that: 

(1) if F K has a unique fixed point Xq then necessarily Xq < min(x_ , xq); 
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(2) since F' K (x) = -jA^J^x) necessarily F' K (x) < 1 for x < x_ ; 

(3) the critical value A c (k) of A such that F K starts to have three fixed 
points is increasing in k. In particular, 

b b 



A c («) > A c (0) = Acr(l) 



(6 - l) fe + x 



(4) if F K has three fixed points then necessarily x_ < x_ and xq < 



(«) / («) 
^ + ' 

(5) the smallest fixed point x_ is continuously differentiable in k > 0. 

Indeed, by the implicit function theorem and the fact that F K {x_ ) < 
1, it follows that 



dK £F K (x M ) - 

_J_ F /■(«)> 






(«) 



A(l-J^(a; (K) )) + K' 
In what follows, let 

(12) m := liminf X n and M:=limsupX„. 

n n 

We are now ready to prove our first result. 

Proposition 3.1. Assume that A > is such that F\ has a unique fixed 
point. Then M = m and hence the recursion (8) has a unique fixed point. 

Proof. Since Y n > 1, it follows from (9) and (11) that X ra +i < Fi(X n -i). 
Since F\ is S'=shaped and is assumed to have a unique fixed point, this 
implies that M < Xg . Moreover, recalling that m = liminf n X n , we see 
that for any e > 0, X n > m + e for all n large enough. Hence, (9) and (11) 
imply that for all large enough n, X n+ \ < F K (X n -i) with k = (1 + m + e) b . 
Thus we obtain 

(13) Me(0,x { l +m+e) "). 

Indeed, if F K has a unique fixed point, then (13) follows immediately. On 
the other hand, if F K has three fixed points then we immediately have M € 
(0,x!l) U (xq,x+). But M < Xq and so M < x (by property (1) of the 

-F K 's) and also in this case xq < x (by property (4) of the F K 's), giving 

(13). 

Since e > is arbitrary in (13), we have in fact M < x_ . Similarly, 
using (10) and (11), we see that m > x_ . We want to conclude that 



14 D. GALVIN, F. MARTINELLI, K. RAMANAN AND P. TETALI 

necessarily m = M. We write 



M -m < 



M 



ds 



d 
ds 



,((i+*) 6 ) 



and the sought statement will follow if, for example, 



sup 

m<s<M 



d_ 
ds 



,((l+s) 6 ) 



< 1. 



(«) 



By properties (1) and (4) of F K it follows that x_ < xq, and hence property 
(4) of J2 implies J^(x_ ) < 1. When combined with the expression for 

(k) 

dx_ /dn given in property (5) of F K , this implies that 



d_ 
dn" 



» 



< 



and hence that 



sup 

m<s<M 



i(l+s) b ) 



ds 



M±+s) b ) 



< sup b — < bx 

m<s<M 1 + 5 



(1) 



(k) (1) 

where the last inequality uses the fact that x_ < x_ for any k > 0. 

„(i) 



(14) 



Thus we have to show that x_ < 1/b. For this purpose it is enough to 

bX 



show that Fi(l/b) < 1/b. We compute 

A 



6^1(1/6) 



(1 + A)r 1+ A 



(1+^) 6 



Next, we observe that the map A 1— > 6A/(1 + \ , b ) b achieves its maximum 



at A r 



(i+i) 6 

(1 + l/b) b /(b - 1), where it is equal to -^ 



b 2 -i 



. The latter 



expression is decreasing in b for b > 2 and for & — 2 it is equal to 



16" 



Therefore, if A/(A + 1) < -j|, i.e., A < 8, then the r.h.s. of (14) is strictly 
less than one. We now examine the case A > 8. We write 



bX 



< 



Xb 



be b be b 
< -^-r- < ^^ < 1, 



(1 + (I^) 6 " (1 + 7? ~ A&_1 " 86_1 

Finally the case 6 = 2 and A > 8 is handled directly: 
2A 2A 16 



for b > 3 . 



:i + 



(i+l) 2 ' 



vl + f) 2 



< 



:i + f) 2 



0.77. 



□ 
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Notice that in the proof of the inequality x_ < 1/6 we did not use the 
hypothesis that F\ has only one fixed point. Moreover, we proved something 
slightly stronger, namely 

(15) there exists e(6) > such that 1/6 — x_ > e(6) for any A. 

The following monotonicity property is an immediate consequence of 
Proposition 3.1. Recall that A cr (l) = b b /(b — l) b+1 is the phase transition 
point for the usual (C = 1) hard core model. 

Corollary 3.2. For every A < A cr (l), the C = 2 multi-state hard core model 
has a unique Gibbs measure. 

Proof. If J2 has only one fixed point then the same is true of F±. By Propo- 
sition 3.1 there is then only one fixed point for the recursion (8). The result 
then follows from Lemma 2.2 and Proposition 2.3. □ 

The next result shows that the phase transition for C = 2 is first order. 
Recall the definitions of M and m given in (12) and let e(6) be as in (15). 

Theorem 3.3. If m ^ M then M - m > e(6) > 0. 

Proof. Suppose m ^ M. From Proposition 3.1, it then follows that F\ (and 

a fortiori J 2) has three fixed points x_ < Xq < X, , with x > xq. 
We now show that xq > 1/6. Indeed, since J(xo) = xo and J is strictly 
decreasing, it is enough to check that J(l/b) > 1/6 or, equivalently, that 

V(! + l) b > 1 / b - But A > (h_f)M-i and clearly 

b b b 2b 1 

> 



( 6 _ 1)6+1(1 +1)6 (6 2_l)6( 6 _l) b' 

Since 1/6 — x_ > e(6) by (15), this implies Xq — x_ > e(6). 

Next, since X n = [1 — P n (a r = 2)] _ — 1, we infer that X n is maximized 
by the empty b.c. and minimized by the full b.c. if n is odd (and vice versa if 
n is even). Thus, using the recursive inequality X n+ i < Fi(X„_i), we obtain 
for any odd n, the inequality X n < U n , where {U n ,n odd} is the sequence 
that satisfies the recursion Un+2 = Fi(U n ), with U\ = 0. In particular, 
m < x_ < -r — e(6). If now M < m + e(6) < x Q then necessarily X n < Xg 
for any n large enough and repeated iterations of X n+ \ < iq(X n _i) imply 

M < x_ . At this stage we are back in the framework of the proof of 
Proposition 3.1 and m = M, resulting in a contradiction. □ 



4. The Large 6 Asymptotic Regime 

In this section we set up and then analyze the recursion for any value of 
C when 6 is large. In what follows, e = exp(l). 
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For any j < C set j* = C - j. Also, for A < 1, set A x = £~ A* 
(1 — A)" 1 . Iterating (5) we obtain 



n / •■* \ b 



£ »'E- 



R b Jk) ) 

^n+2u) = ! + — : r 



J / l \k 

E*ME- 



iS fife) 

i=0 \fc=0 nV y , 



In turn, this implies that 



^■• +1 E 



- 1 A* V 



i4 + 2(j) < 1 + V " 

r J _1 \fc C \fc 

v "\ A v "v A 

1+ **•♦' 



< 1 + 



^fc=0 ij6 (fc) 

1 + A" 1 — r— N 



Therefore, by letting X n (j) = R n (j) — 1 we have 



(16) X n+2 (j) < AlX*-t +1 Jp\x n ) = F%\x n (j)) 



where Xj := A x A J , and j( x > = J and j\ — J2 are the maps defined 
(6), but with the A dependence now denoted explicitly. 



m 
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In a similar fashion, we obtain a lower bound 

\ b 



Rn+2(j) > 1 + - 




1 + 



v^C A fc 



E 



J-l A fc 



fc=0 fl*(fc) 



> 1 + 



^ +1 






'R b nU) 



Therefore, we have 



(17) X n+2 (j) > A- 2 \i*^+ l J^\x n ) = F { l\x n {j)) , 



where A' := A\XJ . 



4.1. The case of C odd. We start by stating the main result of the section. 
Recall that for A < 1, A x = (1 - A) -1 . 

Theorem 4.1. Let j c = |~^-~|, and define A_ := A^ l \^ c and A + := A x X jc . 
Then the following two properties hold: 



i 



(1) If (?) jc < A < 1 with 7 > e , £/ien, /or any 6 £arge enough depending 
on 7, £/ie smallest fixed point of 

(18) x^A\4 X -\x) 

is strictly smaller than the largest fixed point of 

(19) x^A- 2 4 X+ \x). 
In particular, there is phase coexistence. 

(2) On the other hand, if A < ( 3r J Jc u>i£/i 7' < e i/ten, /or every b large 
enough, depending on 7', i/iere is a unique Gibbs measure. 

We start by establishing the first assertion of the theorem. Our proof will 
make use of the following elementary observation. 

Lemma 4.2. For 7 > the function H^ : [0, 00) *—?■ [0, 00) defined by 

H 1 {z)= 1 e- ie "\ z€[0,oo), 
is S-shaped. In addition, the following two properties hold: 
(1) if 1 < e £/ien ff 7 /ias one fixed point z < 1; 
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(2) if 7 > e then H^ has three distinct fixed points z_ < z$ < z + that 
satisfy 

(20) < z- < log( 7 ) - log(log(7)) < z < log 7 < z + . 

Proof. The function H^ is clearly twice continuously differentiable, satisfies 
i? 7 (0) = 7 > and sup^ H 7 (x) = je~ y < oo. That it is S'-shaped therefore 
follows from the fact that 

H'^z) = je- z H y (z) > and Hfa) = je^H^z) [je~ z - l] . 

Now suppose 7 < e. Then sup z H'(z) < 1 and therefore there exists a 
unique fixed point zq. The fact that zq < 1 follows from the observation 
that 

jj 7 (l) = 7e -T e_1 <1. 

On the other hand, if 7 = e the value zq = log 7 is the unique fixed point, 
and satisfies H' (zq) = 1. Lastly, for 7 > e, we have the inequalities 

A 7 (Iog7) > 1, 

#7( lo S7) > lo S7, 
^7( lo g7-l°g(log7)) < log7-log(log7), 

where the last inequality holds because if-,, (log 7 — log (log 7)) = 1 and 
7 I—)- log 7 — log (log 7) restricted to the interval [e,oo) is increasing with 
log(e) — log(log(e)) = 1. Together with the 5-shaped property of H, these 
inequalities immediately imply that H has three fixed points that satisfy 
(20). ' ' □ 

We are now ready to establish the first statement of Theorem 4.1. 

Proof of Theorem 4-1(1)- Fix A € [(r ) Jc ,1) with 7 > e , and for notational 
conciseness, denote A\ simply by A. We first show that the asserted inequal- 
ity between the fixed points of the two maps implies phase coexistence. This 
is a simple consequence of the fact that, for any boundary condition r, the 
sequence {X*} defined by 

X* n = X n (j c ) = HT n (Vr > Jc)/VT n (. a r < 3c), n € N, 
obeys the recurrence 

A- 2 4 X+) (X n ) < X* n+2 < A 2 4 X -\x* n ) 

where we have made use of (16) and (17), together with the duality property 
j* + 1 = j c . If now L^-J boundary conditions are imposed at the zeroth level 
then Xq = and X* will always be smaller than the smallest fixed point 

of x i— y A 2 J \ ~ (x). On the other hand, under |~^] boundary conditions, 
Xq = 1 and X* will always be larger than the largest fixed point of x *— > 

A^J^ + '(x) because the range of this mapping is contained in [0, 1] for large 
b. 
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We now prove our statement concerning the fixed points of (18), (19). 
First, consider the case A = (?) jc and observe that for any z > , 

(21) lim bA~ 2 4 X + \z/b) = lim bA 2 J^ ~\z/b) = H^z), 

uniformly on bounded intervals. Next, we define 

log7-log(log7) _ log 7 

:r_ := and x + := . 

b b 

From Lemma 4.2, it follows that H^ibxJ) < bx~ < bx + < H^(bx + ). To- 
gether with (21), this shows that for any b large enough, 

A 2 4 X -\x-) < x_ < x + < A- 2 4 X+ \x+) , 

and the first assertion of the lemma follows (for this case) because A~ 2 J}, 
and A 2 J \ are S'-shaped exactly like H^. 

We now consider the case (x) jc < A < 1 and again we compute 



' - = A 2 X 

-bx- \ b / .— - \ h ' 



(22) A 2 4 x ~\x-) < A 2 , = A 

(l + \-e- b *-) b ( 1 + X _^ 

If A does not tend to zero as b — )■ oo, then it is obvious that the r.h.s of (22) 
is smaller than X- for large enough b. If instead lim^oo A = we proceed 

as follows. The function / 7 (A) = A/ ( 1 + A- 2 ^ J satisfies 

1 (^ 6Alog7 

1 A- \l 



J:iX} ' lo g ,\ 2b V -+Alo S - 



7 

and hence is decreasing in the interval ( / fe _ 1 7 lo — , oo)- Since 7 > e and our 
assumption A —¥ implies A = A\ ~ 1 for large b, we have the inequality 

A_ > A'^/b > 7/((6 - 1) log 7). 

Thus, we can conclude that the r.h.s of (22) is smaller than the same ex- 
pression with A_ replaced by A^^-j/b. After this replacement, the resulting 
r.h.s of (22) is indeed smaller than x_ for all large enough b because of (20) 

and (21). In conclusion, we have shown that for any (■^) jc < A < 1 the 
function A 2 4 has a fixed point smaller than x_. 

Next, we examine A~ 2 4 ■ If li m b->oo b\+ = 00 then it easily follows 
that for large b, we have A~ 2 J^ A+) (,4~ 2 A + /2) > A~ 2 \ + /2 > x— If instead 
A + < C/b for some finite constant C, we choose x\ = log(6A + )/6 > x_ and 
write 

A- 2 4 X+ \x x ) > A- 2 \ + e~ bX +K l+x ^ b . 

By construction, lim^oo e~ +' ( 1 + x x) _ e -i_ Therefore, for sufficiently 
large b, 

A~ 2 \ + e ~bx + /(i+*x) b > (1 _ 0(b~ 1 ))X + e- 1 > x x , 
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because A+ > 7/6 with 7 > e. In conclusion A~ 2 J 2 (x) has a fixed point 
strictly bigger than x_ and the existence of a phase transition follows. □ 



We now turn to the proof of the second assertion of Theorem 4.1, namely 

the absence of a phase transition for A < I y 1 3c , with 7' < e. For this, 
we first establish two preliminary results in Lemmas 4.3 and 4.4. For any 
vertex y G T n and i G Sp, we define a probability measure fly on the set of 
spins at y as follows: 

(23) ^)(a y = j)=¥(a y = j\a y <i*), j G S c , 

with P, as always, depending on A and a boundary condition on T n (which 
for clarity we have suppressed in the notation) . Note that if a; is a site in T n 

that is neighbouring to y, then fx y represents the marginal on y of the Gibbs 
measure (with some boundary condition on the leaves of T n ), conditioned to 
have i particles at x. Recall that || • ||tv denotes the total variation distance. 

Lemma 4.3. For any k < i, we have 



H'v H'v \\TV 



4'KeP + MD 



lA) (°y < k*) 

Proof. By definition /j, y (a y = j) = fi y (o~ y = j\o~ y < i*). Therefore, also 
recalling that k < i implies k* > i* , we have 



2 Z) II $ ( a v =■?')-# K = i) || + 5 J2 4 fe) ( a y = o) 

j=0 j=i*+l 

l /4 0) K<fc*)-My 0) K<»*) l 4 0) (i* + l<<T y <A;*) 

2 4 0) K<^) 2 /4 0) K<fc*) 



A*i (o"j/ < ^* 



□ 



Notice that if x is an ancestor of y then //j, is nothing but the Gibbs 
measure on the tree T^ rooted at y with the boundary conditions induced 

by those on T n . If instead y is an ancestor of x then fj, y becomes a Gibbs 
measure on the (non regular) tree T n \ T b x . However, if x, y are sufficiently 
below the root of T n , then T n \1* will coincide with a regular tree rooted at 
y for a large number of levels. That is all that we need to prove uniqueness 

below (|)^. 

In what follows, given any non negative function b 1— >■ f(b) of the degree 
of the tree T 6 , we will write f(b) ~ if linib^oo bf(b) = 0. 
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j_ 

Lemma 4.4. Fix 7' < e and assume A < ( -5- ) 3c • Then there exists a < 1 
and no € N such that for any n > hq and any boundary condition r on the 
leaves ofT n , 



6^00 a 1/ Z = J c = I f j 



limsup6 / [/ r (<7 r > i* + 1) < 



Proof. It suffices to bound X n (i) from above for i < |_^-J or i = |~^-~|. In the 
first case, when i < [-j J, the stated bound follows easily since (16) and the 
assumed bound on A imply that for some finite constant K, 

bX n {i) < X l * +l b < Kb^"^ wO. 

In the second case, when i = j c , set Oqo := limsup^^ bx+(b), where x + (b) 
is the largest fixed point of the S'-shaped function x 1— > A\J^ ~ (x). Due to 

the assumption A < ( ^-J 3c , it follows that Oqo < 7'. Because of (16) it is 
enough to prove that aoo < 1. Assume the contrary. Then the fixed point 

equation, together with A < ( %- ) 3c , readily implies that 



< 7e" 



-'ye 



which in turn implies that a^ must be smaller than the unique fixed point 
zo of the map H . Since r ye~' y ' e < 1 if 7 < e necessarily zq < 1 and we get a 
contradiction. Note that in the above proof by contradiction, the hypothesis 
ffloo > 1 enters as follows. Ifx> 1 — 5, < 5 <^l I then Jr, (x) is increasing in 



1 

_Y 

6 



A and so we may safely assume A = ( %- ) 3c and not just smaller or equal. □ 



We are now ready to prove uniqueness for A < 



b 



1 



Proof of Theorem 4-1(2). For simplicity we begin with A = I %- I Jc . In this 

case, it follows immediately from the basic inequality (16) that for any initial 
condition, any n > 2 and any b large enough, there exist constants ci, C2 
such that 

(24) 4 C/2l <c ie - C2ba , 

where a = l/(j' c + 1). In another words, recalling the probability measure 
Hy introduced in (23) and using the obvious fact that for any i < C , 

^(\j c + l,C])<Xl c / 2 \ 

we get that the probability of having more than j c particles at y given i 
particles at x is exponentially small in b. 

Now, recall that 7) is the finite-tree of depth I rooted at r, and let r, r' 
be two boundary conditions on the leaves of T^ that differ at only one vertex 
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vo. Let also T = {vo,v±, . . . ,vg} be the unique path joining vq to the root 
r = V£. We recursively couple the corresponding measures [i T = /jJj, A and 

fi T = fjJp A by repeatedly applying the following step. Assume that, for 

any pair (a Vl ,a V2 ) with a vi / a' vi we can couple /x r (- | a Vl ) and fi T (• | a' Vl ) 

and call Vi-y V1 the coupled measure. It is understood that Ve_y n is 

concentrated along the diagonal if a Vl = a' Vl . Let ir-^ ' v ° be the coupling of 
the marginals on of the two Gibbs measures on v\ that realizes the variation 

distance (i.e., 71"/°' v °(<r Vl ^ a' Vl ) = ||//J — MuJItv)- Then we set 

<Tv l' <T v 1 I l\ Tv 0> T v , / \ <Tv l' <T v 



[a, a = 7r 






If we iterate the above formula we finally get a coupling v T,T such that the 
probability of seeing a discrepancy at the root can be expressed as 

(25) £ vr[- ' T - K , < ) C 1 A fa* , v' V2 )*T '^ • • • 

Vv2^Vv 2 

with self explanatory notation. If we can show that the above expression 
tends to zero as £ — )■ oo faster than b~^ uniformly in t,t', then uniqueness 
will follow by a standard path coupling (or triangle inequality) argument 
(see, for example, [6]). 

On the state space S := [0, . . . , C] 2 consider a non- homogeneous Markov 
chain {£t}t=o with transition matrix at time t given by i"t(£, £') = tT((£') 
and initial condition £o = ( T v , t' v ) . Let also B = {(i,j) € S : i > j c + 
1} U {(i,j) e S 2 : j > j c + 1} be the bad set and let D = {(i,i) G S : 
i € [0, ...,C]} be the diagonal. Equation (25) is then nothing but the 
probability that the chain does not hit D within time £. 

For b large enough (depending only on 7' < e) the two key properties 
of the chain, which immediately follow from Lemmas 4.3 and 4.4 and the 
inequality (24), are the following: 

(26) sup sup P t (tD c ) < £, a < 1 

t £eB c 

(27) sup sup P t (£,5) <cie- C2ba , a >0 . 

Notice that it is not difficult to show that 

sup sup P t (£,D c ) « A»l/6. 

In other words, the probability of not entering the diagonal D in one step is 
suitably small (i.e., smaller than a/b, a < 1) only if we start from the good 
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set B c . Using (26) and (27), we can immediately conclude that 



c (& i D for all < t < £) < ^T ( J 



fc=0 



28 < - -cie- C2b +- 

a \a o 

The "—1" in the exponent of a/b above takes into account the fact that we 
may start at xq in the bad set B, while the extra "— k" in the exponent 
accounts for the fact that for any transition from B to B c we do not nec- 
essarily have a good coupling bound. It is clear that the right hand side of 
(28) tends to zero faster than b~ e as £ — )■ oo because a < 1. □ 



4.2. The case of C even. Throughout this discussion, we assume C even 
§ + 1. Notice that j c = {% 



and we set j c = % + 1- Notice that j c = (§)* + 1. 



i 



Theorem 4.5. ylssitme A = (7-^f- ) Jc wi/i 7 > 1/(C + 2). Then, for any 

large enough b there is phase coexistence. If instead 7 < -r^ro, for any large 
enough b there is a unique Gibbs measure. 



Proof. Fix 7 > ^tt^ and assume A = I 7— P- J c . We will show that the 

largest fixed point of F_ 2 is strictly larger than the smallest fixed point of 

(-) 
F, 2 . By the usual argument that is enough to prove phase coexistence. 

(-) 1 b 
Pick a halfway between l/(C+2) and 7 and compute the value F_ 2 ( Q ? g ) 

for large b. From the definition we get 

(§) alogfe 7 lQ gfr -6^T2- Q 7 log b a log b 
~ l b >~ b C ~ b b ' 

(-) 1 b 

Therefore there exists a fixed point of FJ 2 greater than Q ? g . On the other 

hand 

(29) F® (2^e- bmC+2) ) « 7^e- fe ^ « 2 7 ^e- bl/(C+2) , 
\ J b b 

so that -F + 2 has a fixed point smaller than 27-2|-e , and now the 

first statement of the theorem follows. 

Assume now 7 < ttW. In that case, using (16), we infer that, for any 
boundary condition and any large enough b, 



»T n K > % + 1) < 4 C/2) < e" 



6" 



C + 2 ' 

The proof of uniqueness follows now exactly the same lines of the odd case 
with the difference that now the bad set is B = {C/2 + 1, . . . , C} and (26), 
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(27) are changed into 



(30) 
(31) 



-c 2 b° 



sup sup Pt(£, D c ) < c\e 
sup sup P t (£,£) < cie- C2b ° 



a > 
a > 0. 



□ 



4.3. First-order phase transitions for C even and large b. We now 

turn to showing that for all even C and large enough b (depending on C), 
the phase transition established in Theorem 4.5 is first-order. At the end of 
Section 4.2 we showed that as A varies, for example, in the interval 

2 2 

log6\c+2 /31oe;6\c+2 



the values of 

m(A) := limsup [fS^r > C/2) - ^(ov > C/2)] 

vary between and fi(2-^£-) . (Recall that the superscripts C and indicate 
full b.c. and empty b.c, respectively.) Notice that, by monotonicity, the 
limsup n above is attained over the sequence of even n's and that ^ 2n (a r > 
C/2) is decreasing in n. 

Here, we argue that in the above interval m(A) cannot be continuous. 
The starting point is the observation that, because of (29), for all 



AG 



102-6X^+2 / 3 w b\c+2 



,(¥) ,. 



the smallest fixed point of F + 2 is exponentially small in b a for some a > 0. 
Thus, in particular, there exist constants c\ , ci such that 

H° Tn (a r > C/2) < Cl e' C2 b " , Vn > 1 . 

Fix now 5 < 1 and assume that for some no, 



,C 



S 



u^ (a r > C/2) < - . 
By monotonicity that implies 

sup sup u£ n (cr r > C/2) < - . 

n>2no t 

Thus we can proceed with the previously described coupling argument with 
(30) and (31) replaced by 

(32) sup sup P t (t,D c ) < j, 

t>2n £eB c 

sup sup P t (£,B) < c ie - C2ba . 

t>2n ^ 



(33) 
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and we may conclude that m(A) = 0. 

In other words we have shown that m(X) > implies that for all n, 

5 



so that 



/4> r > C/2) > b , 



m(X)>^- Cl e- C2ba . 
It follows now that the phase transition is first-order. 
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